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Abstract. Lattice molecule models are proposed in order to study statistical 
mechanics of glass transition in finite dimensions. Molecules in the models are 
represented by hard Wang tiles and their density is controlled by a chemical potential. 
An infinite series of irregular ground states are constructed theoretically. By defining 
a glass order parameter as a collection of the overlap with each ground state, a 
thermodynamic transition to a glass phase is found in a stratified Wang tiles model on 
a cubic lattice. 



PACS numbers: 64.60.De,05.60.+q,75.10.Hk 
1. Introduction 

The equilibrium statistical mechanics successfully describes various types of phase 
transitions including ferromagnetic-paramagnetic transitions, gas-liquid transitions, and 
liquid-solid transitions. For structural glass transitions, whose precise definition is 
not obvious, the understanding has been accumulated from several viewpoints pj. In 
particular, in addition to an insightful phenomenological argument [2], [3l HJ [5] , which is 
often referred to as a random first- order transition scenario (RFOT), the analysis within 
the equilibrium statistical mechanics has provided quantitative results for structural 
glass transitions |6l [71 El El [TOl |1T]. At present, it has become widely conjectured 
that a thermodynamic glass transition, if it exists, is described as a one-step replica 
symmetry breaking (IRSB) in the spin glass terminology. Despite such successes, a 
theory of glass transition in finite-dimensional, short-range interaction systems is still 
one of challenging problems in physics, because the theory on the basis of equilibrium 
statistical mechanics has been established only for models with infinite-range interaction 
or on a random graph. 

Here, let us recall a history of studies on critical phenomena. The van der Waals 
theory is the first breakthrough of theory for critical phenomena, which might correspond 
to RFOT for glass transitions. It should be noted that the statistical mechanics 
of a model with an infinite-range interaction exactly predicts critical phenomena in 
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accordance with the van der Waals theory [12]. This reminds us a relation between 
IRSB in mean-field type models and RFOT in phenomenology. Then, the existence of 
a critical point within the equilibrium statistical mechanics of short-range interaction 
Hamiltonians was shown by Peierls |13], Kramers- Wannier [T^, and Onsager [15]. In 
particular, Onsager solved the two-dimension Ising model exactly and proved that 
critical phenomena in finite dimensions are qualitatively different from the van der Waals 
theory (or the Curie- Weiss theory in ferromagnetic-paramagnetic transitions). Since 
then, the significance of critical phenomena in finite dimensions has been recognized 
and much effort has been done in order to connect between Onsager's result and the 
van der Waals theory. We now understand a great picture of critical phenomena. 

However, with regard to glass transitions, there is no exactly solved example in 
finite dimensions; rather, there is no finite-dimensional, short-range interaction model 
for which the existence of a glass transition is understood theoretically. Thus far, toward 
establishment of statistical mechanics of glass transition in finite dimensions, several 
lattice models have been proposed. One example is a class of models proposed by Biroli 
and Mezard [8], in which at most k neighboring particles are allowed to contact with 
each particle. Although this simple model exhibits a glass transition when it is defined 
on a random graph, crystallization occurs in finite-dimensional lattices except for subtle 
cases [16]. In the other model proposed by Ref. [17], crystallization might be prohibited, 
and the numerical experiment has been performed in order to explore the nature of 
thermodynamic glass transitions [18]. However, it seems difficult to develop a theoretical 
argument for this model in finite dimensions. Furthermore, finite-dimensional quenched- 
disordered spin models that exhibit IRSB under the mean-field approximation have 
been studied numerically [19]. However, the numerical computation is much harder 
than standard spin glass models, and a precise theory for the model might be quite 
challenging. (See Ref. [20] as such a theoretical attempt.) 

In contrast to previous studies, we first propose finite-dimensional hard-constraint 
models for which ground states can be constructed theoretically. Here, the ground 
states are obtained by taking the limit of the chemical potential to be infinitely large, 
because the chemical potential is the only thermodynamic intensive variable in such 
models. Indeed, we can show that our model possesses uncountably infinite number 
of ground states in the infinite size lattice, and typical ground states are irregular in 
the sense that they do not exhibit any long-range positional order characterized by the 
existence of Bragg peaks. Note that quasi-periodic ground states are classified as regular 
ground states. Let us denote a set of all ground-state configurations by V, which we 
can specify theoretically for our models. Now, we introduce an overlap with a ground 
state configuration a & V, which is denoted by qa- Since is defined for each ground 
state a, we have an infinite- dimensional vector q = {qa)aev- We refer this vector to as 
the order parameter, because this corresponds to the magnetization in the Ising model. 

We explain the correspondence by reviewing the phase transition in the two- 
dimensional Ising model. When the temperature is higher than the critical temperature, 
in the thermodynamic limit, there exists a unique expectation value of observables with 
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being independent of boundary conditions. However, the independence of boundary 
conditions is broken below the critical temperature. In general, a state of the system 
without the uniqueness is referred to as the ordered phase. The dependence on boundary 
conditions is most easily observed when the expectation of the magnetization is 
considered under the spin-up boundary condition or the spin-down boundary condition. 
We here notice that the magnetization is equivalent to the overlap with one ground 
state. We thus consider the overlap with each ground state as the generalization of 
the magnetization. We also generalize spin-up and spin-down boundary conditions to 
special boundary conditions that uniquely determine a ground-state configuration for 
a series of system sizes going to the infinity. We call such a boundary condition a 
GS-boundary condition. 

Since we have the glass order parameter q in our models, we can investigate whether 
or not the expectation value of q takes a different value under every GS-boundary 
condition. If the dependence is shown, the existence of an ordered phase characterized 
by the parameter q is claimed. From a fact that a typical ground state does not exhibit 
a long-range positional order, we identify the ordered phase as the glass phase. As far 
as we know, such an approach to glass problems has never been attempted. 

This paper is organized as follows. In section |2l we start with a definition of lattice 
molecule models we study. The molecules in the models are represented by hard Wang 
tiles [21] and the molecule density obeys a grand canonical ensemble. In section [3l we 
consider a simple model in this class. We show that the model possesses an infinite series 
of irregular ground-states, while no thermodynamic transition occur. In section HJ as an 
extension of the model, we propose a three-dimensional model in which a thermodynamic 
glass transition is observed. We present theoretical arguments and numerical evidences 
for the thermodynamic glass transition. Section [5] is devoted to concluding remarks. 



2. Lattice molecule model 



Figure 1. Unit square tile witli colored edges. 



Let Al = {(hj) G N X N|l < i,j < L} be a square lattice. We formulate a 
statistical mechanical model in the lattice A^. Each site can be occupied by at most 
one molecule. A molecule is characterized by its shape, represented by m-colors given 
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on edges of a unit square. Since tiles with colored edges are called Wang tiles, the 
molecules in our model may be interpreted as Wang tiles. 

When left, right, bottom and top edges of a Wang tile are colored by £, r, b and 
t, respectively, we denote the quartet of colors by {i,r,b,t). Below, the colors will be 
identified with integers if the correspondence is explicitly given. (See figure [H) Among 
different Wang tiles, we select p different tiles. These are called prototiles and the set 
of prototiles is denoted by C. Each prototile is represented by an integer n, 1 < n < p. 
At each site x G A^, we define a{x) = n if there is a molecule congruent to a prototile 
n E C and as (j{x) = if the site is empty. We call an empty site a hole. The set 
of tile configurations {C U {0})^^ is denoted by S^. We study statistical mechanics 
of the molecules under a boundary condition imposed at sites in (9A^ and 9A^, where 
dKl = {{i,j) G Alii = 1 or j = 1} and dAl = {{ij) G A^li = L or j = L}. The bulk 
region Ai^\(9A^ U dAj^) is denoted by A^,, and the number of sites in the bulk region is 
N={L-2f. 

Specifically, we consider the case that the interaction between molecules is described 
by a hard constraint that molecules are allowed to contact each other only when 
contiguous edges of tiles have the same color. We then assume a grand canonical 
ensemble 

P,i^) = ^^D{a)e^^^>'^^\ (1) 

where D{a) = 1 for configurations that satisfy the constraint that contiguous edges have 
the same color, otherwise D{a) = 0; is the chemical potential of molecules taking a 
value in [— oo, oo], and p{a) is the density of molecules defined by 

/'(^) = ^E[l-^(^(^)'0)]- (2) 

Here, the temperature is set to be unity and its value is irrelevant for the problem. The 
normalization constant in ([1]) is the partition function, which is explicitly given by 

= J2 D{a)e^''^^^\ (3) 

Configurations realized in the limit /x — )■ oo are ground-states in statistical mechanics. 
A tile configuration without holes, which is referred to as a complete tiling in this paper, 
provides a ground state. 

Thermodynamic properties associated with the density are determined by the 
pressure function p{fi) defined by 

p{fi) = lim ^logH(/i). (4) 

The expectation value of the density (p) is given by 
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Furthermore, the entropy density s is related to the pressure in term of the Legendre 
transformation 

s{p) = mi[p{fi) - pfi]. (6) 

Statistical behavior of the model depends on the choice of a set of p-types of 
molecules C. As the simplest example, let us consider the set C with p = 2, in which 
all the edges of one type are red and all the edges of the other type are green. In 
this model, ground states are understood as the complete tilings occupied by one color 
when the open boundary condition a{x) = for x G dAj^ U dA~^ is assumed. When 
fi is sufficiently large, the number density of red tiles depends on boundary condition 
even in the limit L — )■ oo. On the other hand, when /i is sufficiently small (negatively 
large), tile configurations are disordered and all statistical quantities are independent of 
boundary conditions in the limit L — )■ cxd. The Z2 symmetry breaking occurs at some n 
beyond which there exists an ordered phase. The universality class near the transition 
is identical to that of the two-dimensional Ising model. 

A unique feature of Wang-tiles is that the operation of any Turing machine is 
simulated by a complete tiling for an appropriate set C. (See Ref. [21] for the research 
history. See also Ref. [22] as an instructive paper for this issue.) According to 
computation theory, this means that there is no algorithm that determines whether 
or not a complete tiling is possible for a given set C [23j. In solving the decidability 
problem, it was a crucial step to find an aperiodic set of prototiles C for which an 
aperiodic complete tiling in Aqo exists, while periodic complete tilings cannot be realized. 
After the first discovery, p, the number of elements of an aperiodic set has been reduced. 
At present, the minimum number of p is 13 [24J. Statistical mechanics of a system 
consisting of an aperiodic set of 16-prototiles was studied in Ref. [25], where holes are not 
considered, but a positive energy is assumed for mismatches of contiguous colors. This 
reference claims that a thermodynamic transition occurs at some finite temperature. 
See also Ref. [26] for the recent study on the model. 

Here, from a viewpoint of statistical mechanics, it is important to recognize that 
there exits a set C with which non-trivial ground states can be obtained as complete 
tilings corresponding to computational processes. We do not need to stick to aperiodic 
sets of Wang tiles. More important thing in the context of glass problems is that ground 
states should not possess any long-range positional order. However, in the construction 
of complete tilings with using the aperiodic sets with p = 13 and p = 14 [27], the 
quasi-periodic maps are employed to yield the tiling [28J, and therefore their complete 
tilings possess the quasi-periodic order. The construction method in the other cases are 
entirely different from the cases that p = 13 and 14, but at least for known examples in 
Ref. [21], the complete tilings seem to exhibit a long-range positional (quasi-periodic) 
order. 

Now, let us recall a dynamical-system theory, which tells us that aperiodic motion 
described as a solution of a deterministic equation is further classified into quasi-periodic 
motion and irregular motion [29| [30] . Periodic and quasi-periodic motion are called 
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regular motion and characterized by the existence of the singular peak in its spectrum. 
Then, there are an infinitely number of periodic orbits in typical chaotic systems and 
the exclusion of periodic orbits can be realized by systems that exhibit quasi-periodic 
motion. Similarly, it is reasonable to classify aperiodic configurations generated by 
a deterministic rule into quasi-periodic and irregular configurations. Quasi-periodic 
configurations are characterized by the existence of a singular peak in the Fourier 
transform (Bragg peak) of some representation of configuration, as is known in quasi- 
crystals [31j. See also Ref. [32] for a mathematical argument of the definition of 
weak crystals which cover certain generalization of quasi-crystals. In our viewpoint, 
non-periodic long-range positional order in Thue-Morse sequences [33] is classified 
into the same group as quasi-periodic order. Here, it should be noted that irregular 
configurations without any Bragg peaks can also be generated by a deterministic rule. 
In order to seek for thermodynamic glass transitions, we study statistical mechanics 
associated with such irregular (no periodic and no quasi-periodic) complete tilings in 

Aoo- 

From this reason, we are not concerned with aperiodic sets of prototiles, but with the 
case that complete tilings are irregular almost surely when a complete tiling is picked 
up with the equal weight, while there exists a countably infinite number of periodic 
complete tilings. We shall provide an example in the next section. 

3. 4-prototiles model 




Figure 2. Set of 4-prototiles. The left figure shows a number representation. The 
tiles satisfy t = r = £ + b (mod 2) . The right figure shows a color representation. 

The simplest set of Wang prototiles that generates irregular complete tilings is 
given in figure O This set is characterized by a rule that the quartet {i, r, b, t) satisfies 
t = r = i + b (mod 2), where r, i, b and t are either or 1. For this model, we shall 
show that (i) there are an uncountably infinite number of irregular complete tilings, but 
(ii) no thermodynamic transition occurs. 
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Table 1. Cellular automaton rule for generating complete tilings in Aqq. When r of 
the left tile and t of the bottom tile are given, the tile is uniquely determined. The 
rule is equivalent to Rule 102 in Ref. [34j- 



3.1. Complete tilings 

Complete tilings in this model are obtained as follows. For an element 
i[h{t)]t„[v{j)]f^,) in the set Vl = {0,1}^ x {0,1}^, we set = v{j) and 

b{i,l) = h{i). Then, by the rule in tabled! tiles at (i,l), 2 < i < L, are determined 
from smaller i in order. Similarly, the tiles at 1 < « < i^, are determined for 

j = 2, 3, ■ ■ ■ , L. In this manner, all the complete tilings in A^, are uniquely coded by Vl. 
That is, there are 2^^ complete tilings for the system of size L. The complete tilings 
in Aqo are obtained in the limit L — )■ oo. Since I^oo has one-to-one correspondence 
with real numbers in the interval [0, 1], the cardinality of the complete tilings in A^o are 
uncountably infinite. Since periodic tilings are countable, almost all complete tilings 
are aperiodic. 

A remarkable feature of the complete tilings is the additivity. Suppose that ai 
and (T2 are different complete tilings, respectively. Here, we define the addition of two 
configurations ai and a2 by determining i{x) and b{x) as 

£{x) = £i{x) + £2{x) (mod 2), (7) 
b{x) = bi{x) + b2{x) (mod 2) (8) 

at each site x G A/,. We denote this addition by o"i ©(72. Then, from the coloring rule in 
figure [2], we find that the configuration ai ©o"2 is also another complete tiling. Note that 
the configuration cr{x) = 1 for all x G A^, is one of complete tilings. This tiling, which 
is denoted by (To,o, is the unit element in the additive group consisting of all complete 
tilings. As the simplest complete tiling other than cxo^o, we generate a configuration from 
a condition that h{io) = 1 and h{i) = for z 7^ io and v{j) = for any j. As displayed 
in the left-side of figure [3l a fractal pattern, which is known as a Sierpinski gasket, is 
obtained. We denoted it by o'i^^o. Similarly, we can generate a fractal configuration aojo 
from f (jo) = 1 and v{j) = for j 7^ jo and h{i) = for any i. All complete tilings in 
Aoo are given by a superposition of these basic configurations as 



[h{i)ai,o ® v{j)ao,j], (9) 



i,n=l 



where £ represents the summation of the configurations in the sense of ©. When we 
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uniformly choose one complete tiling, it is given by a random superposition of Sierpinski 
gaskets, as displayed in the right side of figure |3l The randomness of h{i) and v{j) in 
([9]) yields irregular configurations. 




Figure 3. Examples of complete tilings. L = 32. Left: h{i) — d{i,l) and v{j) — 0. 
Right: h{i) and v{j) are chosen randomly. 

As a preliminary for later arguments, we define the overlap ratio between two 
complete tilings a, a' G V^o- 

qa,a' = J2 ^i(^a{x),aa{x)), (10) 

and we consider the probabihty density Pcsil) that g^.a' takes a value q when a and 
a' are sampled uniformly. Since the 4-prototiles in our model are distributed with an 
equal ratio for a typical sample uniformly chosen from Pqo, we have Pgs(q') = ^il — 1/4) 
in the limit L — )■ oo. 

At the end of section 13. let us remark previous studies related to the 4-prototiles 
model. The complete tilings of the model are essentially the same as the ground- 
state configurations of a three-body interaction spin model defined on upward triangles 
[35| [36l [37] . While an extremely slow dynamics is observed, no thermodynamic phase 
transition occur at a finite temperature. Recently, it has been shown that the system 
under a magnetic field, which is equivalent to a lattice gas model with a chemical 
potential, exhibits a first-order transition in the (T, fi) space [38] . Although non-trivial 
nature of the coexistence phase has been conjectured, further studies remain to be done. 

3.2. Statistical Mechanics 

We consider statistical mechanics of the 4-prototiles model. It is convenient to introduce 
a hole variable 

r]{x) = 6{a{x),0). (11) 

That is, ri{x) = 1 only if the site x is empty. For a given hole configuration r] G {0, 1}^, 
we define a region N^j = {x & AL\ri{x) = 0} in which tiles occupy the sites. Let 
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a : ^ C he a restriction of a on the region N^. Any configuration a can be then 
expressed by variables r] and a. The partition function is expressed as 

H(/i) = e'^^ e-''^-^'^L^^^^Q{r]), (12) 

»;G{0,1}^ 

where ^1(7]) is the number of possible tile configurations a when a hole configuration r) 
is given. We here fix one complete tiling a in Aqo- We assume a boundary condition 
that a{x) = (Taix) for x G and a{x) = for x G A^. We can calculate fl{r]) 
most easily under the boundary condition. Note that the pressure function p{fi) in the 
thermodynamic limit does not depend on the choice of boundary conditions. 

The number fl{r]) is estimated as follows. It is obvious that = 1 when there 
are no holes. Suppose that there is a hole at the site (zi, ni). Then, two tiles are possible 
at the site {ii + 1, ni) if there is no hole at the site (ii + 1, ni). We choose one tile of the 
two. For each choice, another tile configuration is uniquely determined by repeating the 
rule in Tabled] with starting from the site {ii + l,ni). Note that the configurations are 
consistent with the boundary conditions at i = L and n = L. Similarly, there are two 
tile configurations that originate from the tile replaced at the site [ii, ui + 1) if there is 
no hole at the site (2i,ni + 1). Therefore, the number of possible configurations is 

where we set ri{L,n) = ri{i,L) = 1 in accordance with the boundary conditions. 

Next, we assume that there is another hole at the site {i2,n2). Then, in a manner 
similar to the first case, the number of possible tile configurations is assigned to this hole. 
The important thing here is that this estimation is obtained with being independent 
of configurations generated from the first hole at the site (zi,ni). Since there is no 
over-counting of configurations, the number of possible configurations is 

2l-»?(«i+l,ni) . 2l-''?(«i."i+l) . 2l~''(*2+l,n2) . 2l~'?{*2."2+l) (^4^) 

Repeating these considerations for all holes, we arrive at the expression 

By introducing a spin variable s{x) = 2t]{x) — 1, we rewrite ( ITSl) as 

fi(ry) = e^^-^(^), (16) 
with a Hamiltonian 

^(^) = I E ^(^)^(^'), (17) 

{x,x') 

where J = log(2) and (x, x') represents a nearest neighbor pair x and x'. Therefore, the 
partition function given in f[T^ is expressed as 

H(^) = e^^+^^ J2 e-^W-^S^^(^). (18) 

se{-i,i}'^ 
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That is, is determined from the canonical partition function for the anti- 

ferromagnetic Ising model under a magnetic field. Since J/4 = 0.173.. is less than 
the critical point /3c — 0.44... of the Ising model, the pressure function does not 
exhibit any singularities as a function of /i. Thus, there is no thermodynamic transition 
in this model. Related with the problem, we remark that the expression ( ITSll with (|T711 
suggests the existence of a phase transition if J were larger than 4/3c. Such a case 
may arise in a p-prototiles model with sufficiently large p. The transition in this case 
is regarded as an entropy-driven crystallization of holes. Although it is an interesting 
phenomenon, we do not discuss it in this paper. 

Let us recall that there are an infinite series of complete tilings. When /i — )■ cxd 
is considered for fixed L, the expectation value of some observables depends on the 
boundary conditions in the infinite size limit. Note that this limit is different from 
the case /i — )■ oo in the infinite system in which L — )■ oo is considered for the system 
with /i < oo. Since there is no thermodynamic transition, no boundary condition 
dependence of quantities is observed in the latter case. That is, once holes are generated 
in a complete tiling with any positive ratio, the system becomes free from boundary 
conditions despite the infinite degeneracy of complete tilings. We may say that complete 
tilings in the infinite system are unstable with respect to generation of holes. The origin 
of the instability is understood from the fact that an infinite region is infiuenced as the 
result of iteration of the cellular-automaton rule from one hole. (Recall the argument 
above f|T3|) .) Although a cellular- automaton rule can easily generate an infinite series of 
irregular complete tilings, it simultaneously leads to the instability of complete tilings 
so that a thermodynamic transition is not observed. In order to have stable complete 
tilings against generation of holes, we need to avoid a chain of tile-replacements induced 
by one hole. 

4. Stratified Wang tiles 

Let Al,m = j, A;) e N X N X N|l < i, j < L,l < k < M} he a cubic lattice. As 
a natural extension of the models in the previous section, we consider stratified Wang 
tiles in A^^m- In addition to the color matching condition in each plane with k 
fixed, we further impose a constraint condition that two neighboring prototiles in the k 
direction, if they exists, are the same type. We define a{x) = n if there is a molecule 
congruent to a prototile n G C on the site x G A^ a/ and as a{x) = if the site x is 
empty. Specifically, we consider the 4-prototype model studied in the previous section. 
A complete tiling in this model is given by <j{i,j, k) = <7a{i,j) for a G Poo- We denote 
it by aaix) for x G Ax,,m- 

We study statistical mechanics of the model with a boundary condition imposed at 
sites in dAlj^j, dAlj^^, and dA^^lj. where dAlj^j = {{ij,k) G AlmIi = 1 or j = 1}, 
9^L,M = {{^J^ ^ J^L,M\i = L or j = L}, and <9AX% = {{i,j, k) G AL,Ai\k = 1 ot k = 
M}. The number of sites in the bulk region is = (L — 2)^(M — 2). Without confusion, 
the notations in the two-dimensional case are also employed in the stratified model. We 
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assume that a obeys the grand canonical ensemble 

P,ia) = ^^D{a)e^''''^'^\ (19) 

where D{(y) = 1 for configurations that satisfy the constraints that contiguous edges of 
tiles in the {i,j) plane for each k have the same color and that neighboring prototypes 
in the k direction are the same type, otherwise D[a) = 0; is the chemical potential of 
molecules, and p{a) is the density of molecules defined by 

p(a) = l Yl [i-'5K^),o)]. (20) 

The normalization constant S(/i) in ( !T9|) is the partition function of the stratified model. 



4.1. Glass transition 

We explain the existence of a thermodynamic transition in this model. We fix a complete 
tiling a G Poo- We study statistical mechanics (1191) of the system by assuming a 
boundary condition that 

a{x) = <Ja{x) (21) 

for X G dA^ and cr(x) = for the other boundary sites. We refer it to as a GS- 
houndary condition, because the ground-state configuration is uniquely determined by 
this condition. The expectation value under the GS-boundary condition is denoted 
by ( Here, let us recall that the up-spin boundary condition of the Ising model 

under which the ground state is uniquely determined despite the system possesses the 
up-down symmetry. The order parameter of the Ising model is the magnetization and 
it is interpreted as the overlap with the ground-state configuration. Similarly, in our 
model, as an order parameter associated with the complete tiling a, we define an overlap 
variable with a as 

= ^ E (22) 

Now, we consider the case that 1 <^ e^ <^ L = M. A typical configuration in this 
case is given by a random deposition of holes with a probability with + 0(e~^^) for 
each site. When a hole is inserted to the complete tiling, no other tile configurations 
are allowed, which is different from the two-dimensional case discussed in the previous 
section. A tile different from at a site Xq appears when four holes enclose the site 
Xq. The probability of such a configuration is 2e~^'^ -|- 0{e~^^). In general, we expect 
that \miL^oo{qa)a^ is expressed as a convergent expansion in for sufficiently large 
/i. That is, 

lim(g„)^'^ = l-e-'^ + 0(e-2^). (23) 

L— >oo 

for large /x. A more precise estimation on the basis of a Peierls argument should be 
presented. (See a brief discussion in section O) 
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Next, we assume the open boundary condition that a{x) = for all the boundary 
sites. We denote the expectation value under the boundary condition by ( )q'^. Then, 
when /i is sufficiently large, a typical configuration is close a complete tiling, but it is 
not identical to g^, in general. As discussed in section 13. H the overlap ratio between 
two different complete tilings is 1/4. We thus obtain 

hm (g„)o^'^ ^ \ (24) 

for sufficiently large /x. By comparing (123!) and fl24|) . we conclude that the expectation 
value of the observable qa depends on boundary conditions in the infinite size limit. 
This means that there exists the ordered phase in the system with sufficiently large /x. 
In the other limit where /i is sufficiently small (negatively large), statistical properties 
are independent of boundary conditions, because dilute tiles are almost non-interacting. 
Thus, there exists a thermodynamic transition at some value of /i. 

The quantity characterizing the ordered phase is a collection of qa, which is denoted 
by 9 = (?a)aei>cx>- Formally, the order parameter q is an uncountably infinite dimensional 
vector. Such an order parameter is quite peculiar. To our best knowledge, no examples 
in finite dimensions have been reported. We also recall that a typical complete tiling 
a does not exhibit Bragg peaks in its Fourier transform. From these conditions, we 
identify the thermodynamic transition to the ordered phase with the glass transition to 
the glass phase characterized by q. 



4-2. Numerical experiments 

We report results of numerical experiments for the stratified 4-prototiles model. In order 
to facilitate the equilibration, we employ the exchange Monte Carlo method [39]. We 
prepare K replicas of the system with /i^ = fi^^^xk/K, k = 1, ■ ■ ■ , K. In this paper, we 
set /imax = 2 and K = 40. We estimate the expectation value (A)^'*^ for an observable 
A{a) by the time average of A{a(t)) during the time interval [t^,2t^] with discarding 
the transient data [0,t„], where the initial condition was assumed to be cr = cTq, for all 
fik- When the time average of A{a{t)) is independent of tw within statistical errors, we 
assume that the result provides the estimation of equilibrium values (A)^'^'^ . 

In the left side of figure HI p = {p)^^'' is plotted as a function of p for the system 
with four different system sizes. In order to clarify the singular nature, the fiuctuation 
intensity defined by x = ^{{.P — pY)a'^'^ also displayed in the right side of figure HI 
It should be noted that the fiuctuation relation x = dp{p)/dp holds. The graphs of x 
indicate the existence of a thermodynamic transition. 

Next, we numerically investigate the existence of the glass phase by measuring the 
order parameter qa. In the left side of figure [5l we display q = {qa)a^ (9a)o which 
shows the boundary condition dependence of the expectation value of q. From the size 
dependence, we expect that the behavior sustains in the infinite size limit. This means 
that the transition is identified with the glass transition. Note that the expectation value 
of the density is independent of boundary conditions in the thermodynamic limit. In 




Figure 4. Statistical quantities of the density. {L,M) — (12,8) (plus), (15,10) 
(cross), (18,12) (asterisk), and (21,14) (box). Error-bars are within the symbols. 
Left: Expectation value of the density. Right: Fluctuation intensity of density x- 




Figure 5. Statistical quantities of the glass order parameter. {L,M) = (12,8) 
(plus), (15,10) (cross), (18,12) (asterisk), and (21,14) (box). Error-bars are within 
the symbols expect for Xq with {L,M) = (21,14). Left: Expectation value {qa)a'^ 
and {qa)o^'^ (circle). Right: Fluctuation intensity of the order parameter Xq- 



order to characterize the singularity associated with the order parameter, we measured 
Xq = N{{qo, — (lY)a'^^- As shown in the right side of figure [5l Xq exhibits a divergent 
behavior at some value of /x. 

Here, we study the nature of the singularity quantitatively. We first attempt to fit 
the fiuctuation intensities x and Xq with power-law functions as |yU — yud, where /ic is 
a transition point. As displayed in the left side of figure Ej x "^^Y exhibit a power-law 
behavior 

X^I/i-ZicI"", (25) 

with a ^ 0.6 and /ic = 1.175. This result indicates that the transition is the second 
order according to the Ehrenfest classification. On the other hand, a clear power-law 
behavior is not observed in Xq- As shown in the right side of figure El a fitting of the 
power-law form Xq — {f^c ~ fJ')~'^ might be not so bad with 7 ^ 0.8 when we choose 
a value of /ic (yUc = 1-17). However, Xq the regime fi > fic is far from power-law 
behaviors even if we change the value of /ic. The singular nature of order parameter 
fiuctuations is quite unusual. 

In order to demonstrate the peculiar behavior more, we consider the distribution 
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Figure 6. Left: Log-log plot of x as a function of \fic — where fic is assumed to 
be 1.175, for {L,M) = (21,14). The square symbols correspond to the case ji < Hc 
and the circle symbols to ^ > ^c- The guide line represents x — 0-2|(Mc — 
Right: Log-log plot of Xq as a function of |/ic — where /ic is assumed to be 1.17, for 
{L,M) — (21, 14). The square symbols correspond to the case fi < fi^ and the circle 
symbols to /i > /ic. The guide line represents x<? = (a*c — /i)^"'^. 
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Figure 7. Left: Color representation of P^i{q) in the {q, space. {L,M) — (21, 14). 
Right: P^{q) with fj, = 1.15 (square), fi = 1.2 (circle), and 1.25 (asterisk). The error- 
bars are within the symbols. 



function of q for the system with /i, which is denoted by Pn{q)- The left side of figure 
[7] shows that the peak of P^{q) for fi > fi^ accompanies a broad tail in the smaller q 
region. We investigate Pfj.{q) near a transition point carefully in the right side of figure 
[71 Let us focus on the graph with /i = 1.25 represented by the asterisk symbols, which 
corresponds to the second left of circle symbols in the right of figure [6l There is a small 
slope regime around q = 0.5 which is apart from the main peak q = 0.7. The left edges 
of the regime seems to be connected to the left edge of the nearly fiat regime of the graph 
with /i = 1.2, which corresponds to the most left of circle symbols in the right of figure 
[6l The left edge of the flat regime may be further connected to the main peak of the 
graph with fi = 1.15, which corresponds to the most left of square symbols in the right 
of figure [6l On the other hand, the main peak of the graph with = 1.25 arises from the 
right edge of the fiat regime of the graph with fi = 1.2. The important thing here is that 
the main peak in the ordered phase is not connected to the main peak in the disordered 
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phase. This suggests the discontinuous transition of (g)^'*^ in the thermodynamic hmit 
L,M ^ oo, although a clear discontinuous jump is not observed in the left side of figure 
|5l Extensive numerical studies are necessary to have firm evidences for supporting the 
conjecture, because the sizes we investigated are too small. In any cases, we can say 
that the behavior of the order parameter q is quite unusual. 

5. Concluding remarks 

In this paper, we have proposed lattice molecule models, one of which exhibits the 
thermodynamic glass transition in three dimensions. The glass phase is characterized 
by the uncountably-infinite dimensional order parameter q, each component of which 
measures the overlap with an irregular ground state. We conjecture that the transition 
is the second-order in the sense of thermodynamics, while the glass order parameter 
exhibits a unusual behavior. Before ending the paper, we describe problems that should 
be studied seriously in future. 

Although we have studied the system with the simplest set of 4-prototiles that 
generates irregular complete tilings, statistical behavior of lattice molecule models 
depends on the choice of prototiles. By studying models with other sets of prototiles, we 
wish to classify their phenomena. In particular, it is interesting to find a two-dimensional 
model that exhibits a glass transition or to prove that there is no such a model. Note 
that the thermodynamic transition to the quasi-periodic ordered phase was observed in 
a two-dimensional Wang-tiles model [251 126] . The difference between the quasi-periodic 
order and the glass order should be clarified. In doing these studies, extensive numerical 
experiments are necessary. It is significant to develop an efficient method for numerical 
calculation. The hard nature of molecules would substantially reduce the computation 
time if an elegant algorithm is found. 

Since the theoretical arguments reported in this paper are still in the early stage 
of the study, important theoretical problems remain to be solved. The first problem 
is to provide a mathematical proof of the existence of the glass transition. This might 
be solved as follows. We consider a probability that a at the center site xq is different 
from (Tq under the boundary condition a = aa for all boundary sites. We estimate 
an upper-bound of the probability by noting an interface separating an ordered region 
connected to the boundary sites. (See an example of an interface in figure 0) If this is 
explicitly defined for a given arbitrary configuration in which cr(xo) 7^ crai^o), we may 
formulate a Peierls argument. That is, a transition is understood from the competition 
between the entropy cost of configurations with interfaces and the energy cost of holes. 
However, up to the present, we do not have an explicit definition of such interfaces. 

A more important but difficult problem is to obtain a mathematical description of 
statistical properties of the glass order parameter q. When we consider this problem, it 
seems better to forget the tile model. Instead, we will analyze a stratified model of the 
3-body spin interaction on upward triangles. (See a remark at the end of section ISTTl ) 
By applying the same argument in the present paper to the spin model, we may find 
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Figure 8. Schematic figure of a simple interface configuration witli large 77. Holes in 
the top and bottom triangles form Sierpinski gaskets and the vertexes of the triangles 
are connected by holes in the vertical directions. 

a glass transition for the stratified model. Since this model is much simpler than the 
tile model, several theoretical calculations including a sort of the Bethe approximation 
will be done more easily. Such theoretical study may provide a connection to the mean- 
field picture of the glass transition (RFOT). Furthermore, it is amazing if we find an 
exactly solved model that exhibits a glass transition. An exact solved model would play 
a crucial influence on the study of glass transitions. 

In the analysis of infinite-range interaction models that exhibit a glass transition, 
the number of pure states in the glass phase has been one of concerns. With regard to this 
problem, we briefly review Gibbs measures of infinite-size systems. Roughly speaking, 
a Gibbs measure is defined in such a way that the probability of configurations in any 
finite-size region is given by the grand canonical ensemble with boundary conditions 
which are chosen by the measure. In the disordered phase, the Gibbs measure is unique, 
while there are an infinite number of Gibbs measures if the uniqueness is broken. In 
particular, a special measure that cannot be decomposed further into a superposition of 
other Gibbs measures is called a pure state. Here, let us recall that each GS-boundary 
condition in the Ising model can provide a pure state. Although GS-boundary conditions 
are not always related to pure states, the statistical ensembles in finite systems under 
GS-boundary conditions may be a starting point for understanding of pure states in the 
glass phase. 

Finally, we go back to our motivation of understanding the nature of glassy 
materials. Although we have found a thermodynamic glass transition in a short- 
range interaction model in finite dimensions, it is not obvious whether or not such an 
idealized phase is actually observed in laboratory experiments. Toward an experimental 
realization of thermodynamic glass phases, we need to consider the following problems. 
The first is to construct a mechanical model that exhibits a glass transition. Although 
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we have only to design a potential function sharing common features with hard- 
constraint conditions in Wang tiles, its explicit demonstration may be challenging. 
The second problem is to find an experimentally realizable algorithm for facilitating 
the equilibration, because the exchange MC method cannot be employed in laboratory 
experiments. After solving the problems, we hope that we will be able to demonstrate 
by numerical experiments that a genuine thermodynamic glass transition is observed in 
laboratories. 
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